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The so-called generalized Wigner distribution has earlier been shown to be an excellent approx- 
imation for the terrace width distribution (TWD) of vicinal surfaces characterized by step-step 
interactions that are perpendicular to the average step direction and fall off as the inverse square of 
the step spacing. In this paper we show that the generalized Wigner distribution can be derived from 
a plausible, phenomenological model in which two steps interact with each other directly and with 
other steps through a position-dependent pressure. We also discuss generalizations to more general 
step-step interactions and show that the predictions are in good agreement with TWDs derived 
from numerical transfer-matrix calculations and Monte Carlo simulations. This phenomenological 
approach allows the step-step interaction to be extracted from experimental TWDs. 
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X3 

o 



> 

oo 
o 

00 
O 
O 

o 



I 

o 
o 



X 



I. INTRODUCTION 

Vicinal surfaces consist of terraces divided by steps 
which interact with each other via a variety of mecha- 
nisms, including elastic, dipolar, and indirect electronic 
interactions. Since these interactions directly determine 
the distribution of terrace widths, the use of terrace 
width distributions (TWDs) to determine the step-step 
interactions is an important goal. Theoretical attempts 
towards this goal have produced differing results depend- 
ing upon specific approximations P, Q, 0, S IE IE 13 • 

In many cases, the potential energy V due to the inter- 
action between neighboring steps can be written in terms 
of the distance L between the steps as 



(1) 



This model of interacting steps can be mapped directly 
onto the Calogero-Sutherland model |E 13 of interacting 
spinless fermions (or hard-core bosons) in one spatial di- 
mension (Fig. in this picture, the direction along the 
steps (the y-direction in "Maryland notation") is inter- 
preted as time, and the steps themselves are interpreted 
as the world-lines of the spinless fermions. The static 
properties of this system (such as the TWD) depend on 
A only through the dimensionless interaction strength 



A = 



fiA 



(2) 



where (3 is the step stiffness, /cb is Boltzmann's constant, 
and T is the absolute temperature. The TWD can be 
calculated exactly ||,ta] when i = -l/4, 0, or 2. 



'Corresponding author: "Howard'RichardsOtamu-commerce.edu" 

ftp : //faculty . tamu- commerce . edu/HRichards/index . html 
instein@u md.edui, jhttp : / /wwg2 . phy s i c s . umd .edu/'einstein/ 




position 



FIG. 1: Steps on a vicinal surface can be mapped to the 
worldlines of spinless fermions or hard-core bosons in one spa- 
tial dimension (the a;-direction) . The "time" axis corresponds 
to the y-axis on the vicinal surface. 



An argument [lOj from random matrix theory orig- 
inally due to Wigner suggests that these three exact 
TWDs can be approximated by 



P{s) 



gS^ exp 



(3) 



where s = L/{L), the constraint (s) = 1 gives bg [see 
Eq. H32I) ]. flg is a normalization constant, and the rela- 
tionship between A and p is |E 



A 



(4) 



These "Wigner distributions" have proven quite success- 
ful and have been widely used [Til IT^ flSj in nuclear 
physics (in which case s is the spacing between energy 
levels), random matrix theory, quantum transport, and 
universal conductance fluctuations. There have been a 
number of attempts HI EE EE IH E El 

to inter- 
polate between these three special cases. 
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The so-called generalized Wigner distribution j2l| is 
defined by the simple assumption that Eq. J^J) and Eq. Q 
remain valid for all values of A. (This is only reasonable 
for positive or weakly negative values of A The 
generalized Wigner distribution appears to be in bet- 
ter a gree ment with computer simulations of vicinal sur- 
faces [23 than its competitors, and is in reasonable agree- 
ment with many experimental TWDs !25|. Furthermore, 
it has been proved that for rational values of g, P(s) cxs^ 
as s — > [23, l28l| . However, for arbitrary values of 
A there is no clear relation between random matrix the- 
ory and the Calogero-Sutherland model. This paper pro- 
vides a formal and physically intuitive justification for 
Eq. ^ for arbitrary values of A] the resulting formalism 
allows us to consider more general step-step interactions 
than Eq. Some of these results have been mentioned 
briefly elsewhere but are explained here in detail. 

The organization of this paper is as follows. In Sec. ITU 
we show that the generalized Wigner distribution can be 
derived from a phenomenological model, which reduces, 
under certain circumstances, to the two-particle Calogero 
model ||8J. This leads, by separation of variables, to a 
one-variable Schrodinger equation that determines the 
TWD. We extend this treatment to non-trivial examples 
in Sec. Illll The significant process of extracting step- 
step interactions from experimental TWDs is discussed 
in Sec. IIVI Finally, we conclude in Sec. El with some 
possible extensions of this work. 



II. PHYSICAL EXPLANATION OF THE 
GENERALIZED WIGNER DISTRIBUTION 

A. Deriving a Schrodinger equation from the 
generalized Wigner Distribution 

Given a system of steps with interactions of the form 
specified by Eq. and a TWD given by the correspond- 
ing generalized Wigner distribution [Eq. we can de- 
fine a wave function such that iIjq{s)=P{s): 



Differentiating twice, wc find 



a;,'-s-'-exp(-6es72) . 



(5) 



1(1 ^^)'"'"(^+^^^^ + ^'^ ^"^'^ 



As 2 - {g+l)be + bls^ Vo(s) , 



(6) 



where we have used Eq. |@J. Eq. ® allows us, in 
somewhat the same spirit as the Gruber-MuUins ap- 
proximation to propose the following dimensionless 
Schrodinger equation: 



ds2 



where 



^„{s) = [U{s) + V{s)-Enm, (7) 



V{s) = (8) 



is the explicit step-step interaction potential in dimen- 
sionless form and 



U{s) = bl 



(9) 



is a dimensionless potential due to other steps not explic- 
itly considered. (This idea is explained in the next sub- 
section.) By inspection we see that ipois) is the ground 
state wave function (it has only one antinode), with an 
eigenvalue given by 



Eo = ig+ l)bg . 



(10) 



In comparing different experimental TWDs with a 
range of values of {L) but with presumably the same 
step-step interaction V{L), it is often useful to rewrite 
Eq. 121 in its dimensional form 



(fcBr)2 d2 



■i>n{L)^[U{L) + V{L)~En]^jL), (11) 



/3 dL2 
where V{L) is given by Eq. 



and 



E„ 



(12) 



(13) 



Throughout this paper we will alternate freely and with- 
out comment between dimensional and dimensionless 
representations . 

Eq. lO can be solved for all the eigenfunctions and 
eigenvectors |29j . so one can conveniently consider per- 
turbations from the purely inverse-square interaction 
given by Eq. The eigenfunctions are 



V'n(s) 



s^/^exp{-bgS^)iFi 



g+l 



b„s~ 



(14) 

where n is any nonnegative integer, c„ is a normalization 
constant, and if 1 (a, b, x) is Rummer's confluent hyper- 
geometric function l3(l|. The corresponding eigenvalues 
are 



E„ 



{g+l + An)bg 



(15) 



B. 



Deriving the generalized Wigner Distribution 
from a Schrodinger equation 



An undesirable element of the previous derivation was 
that it implicitly required that one of the steps be held 
fixed, which introduced an artificial (and unnecessary) 
difference between steps. Furthermore, the quadratic na- 
ture of U{s) was not well explained. In this subsection 
we explore these issues. 
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The Calogero-Sutherland 's', model with an infinite 
number of interacting spinless fermions would be an ideal 
model of the vicinal surface, but as mentioned above, it 
is only integrable for three values of A. A useful alterna- 
tive is to consider only two adjacent fermions explicitly 
and model the effects of the other fermions phcnomcno- 
logically through the pressure they exert, yielding the 
Hamiltonian 



'H=-- 



1 ( 92 



92 



2 \dx\ 



(16) 

Here Vi is the pressure exerted on fermion i, which we al- 
low to be position-dependent because the limited correla- 
tion length in the x-direction leads to an effectively finite 
system size |3ll |. Accordingly, we place the V fermions 
immediately to the left of fermion 1 in a "box" with a 
fixed left wall at a; = —V; the right wall is at xi, the po- 
sition of fermion 1. The volume (length) of this box is 
Vi = V+xi. The pressure Vi can be expanded as follows: 



Xi 



dVi 
dxi 
dVi 
dVi 



xi=0 



0{xl) 
^Oixj) 



(17) 
(18) 



= V + xi 

= r -xiivny^ + o{xl), (19) 

where 7^ = 7^1(0) and the compressibility k is given by 



K — — — 



dVi 



(20) 



Vi=V 



Likewise, we place the V fermions immediately to the 
right of fermion 2 in a "box" of volume (length) V2 = 
V — X2, yielding 

r2ix2) ^V + X2{Vk)-^ +0{X22). (21) 

Combining Eqs. (fT?)|l . (fll^ . and (|^ . we obtain 

+ {xI+xI){Vk)-^ + 0{xl-xl). (22) 
Finally, we can make the following change of variables: 

Xi + X2 



S = X2 — Xi > 

This allows us to rewrite the Hamiltonian as follows 
'1 92 



(23) 
(24) 



H 



4 dxL 



-1 „2 



Vis) 



+ [Ps + (Vny's^] + A(yK)''xi^^ + 0{sx){25) 



If terms of quadratic and higher order in Eqs. H19|l and 
(|21|l can be neglected, the Hamiltonian can be separated 



into a part that depends only on s and a part that de- 
pends only on Xcm'- 



where 



and 



with 



He 



Ti. — Ti-c 



1 92 



4axL 



4(Vk) 



Us 



+ V{s) + U{s) 



U{s) =Vs+ {Vnyh^ . 



(26) 



(27) 



(28) 



(29) 



None of the argument so far requires the assumption 
that V{s) is given by Eq. (jSJ. In the next section we 
will examine how well this phenomenological formalism 
works when applied to other forms of V{s). 

If (Vk)^^ ^V, we can neglect the term Vs in U{s). 
The resulting Hamiltonian Ti. can easily be seen to be the 
Hamiltonian of the two- fermion Calogero model Q, 



92 



92 



dxl 



dxl 



^ = - hrr + TTT +V{x2~xi) + {x{+x'2)u:' , (30) 



where we identify from Eq. ^ 

= (Vk)-i . 
More importantly, with the identification 



r(^) 



(31) 



(32) 



which again comes from the constraint (s) = 1, we see 
that 



(7^, -^„)i/;„(s) = 



(33) 



is just Eq. Q and again leads to the generalized Wigner 
distribution. 

Finally, it should be pointed out that the one particle 
"Calogero model," with a Hamiltonian given by 



U = -—r+X^UJ^. 

da;2 



(34) 



is really what in essence is used in the many Gaussian ap- 
proximations 

[1 II II II B II II ■ These approximations, 
however, produce conflicting functional relationships be- 
tween Lu and A. This is hardly surprising, since V does 
not appear explicitly in Eq. H34|) . By using a two-particle 
Calogero model, we are able to state unambiguously the 
relationship between lo and A. 
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FIG. 2: [Color online] TWDs for systems of steps with step- 
step interactions of the form V{s) = A^s'^ . Diamonds, circles, 
and squares mark TWDs from Monte Carlo simulations with 
As = 0.5, 2, and 6, respectively. Solid lines are solutions of 
Eq. ^ for the same values of A3 and with U{s)(xs . 



III. NONTRIVIAL EXTENSIONS 

Here and in the remainder of the text, for our numer- 
ical work we use models that are discrete in both the 
y-direction and the x-direction 0) IE US > in contrast 
to the continuum step mode, which corresponds to the 
description given above. Specifically, in Monte Carlo sim- 
ulations we use the terrace-step-kink (TSK) model, in 
which the kinks can be any integral number of lattice 
units long. Our transfer-matrix calculations use the re- 
stricted TSK model, in which the kinks can only be of 
one lattice unit. In both models there are exact expres- 
sions for the stiffness of an isolated step as a function of 
temperature '3^. The interaction between steps is given 
by a specified function V{L), as with the continuum step 
model. Further details of our numerical work can be 
found in Refs.mil 

Eq. (TJ can be extended to step-step interactions of 
forms other than Eq. (jS)). In this section we will present 
a few arbitrary but physically motivated step-step inter- 
actions and show that the phenomenological method of 
the preceding section still yield TWDs in excellent agree- 
ment with numerical simulations of the TSK model. 

As a first example, we consider interactions of the form 

V{s) = Ass-^ , (35) 

which would be a plausible subdominant term in a Taylor 
expansion of the step-step interaction in s^^. Figure [21 
shows TWDs obtained by Monte Carlo simulations of 
steps interacting with this potential for three different 
values of A3. Since the Wigner distribution worked so 
well with nearest-neighbor interactions of the form given 
by Eq. ©, we again take [/(s)ocs^. The solutions of the 
Schrodinger equation are in excellent agreement with the 
Monte Carlo data. 
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FIG. 3: Comparison of a^g from the generalized Wigner 
distribution (solid curve) and transfer-matrix calculations for 
nearest-neighbor-only (NN) interactions (dashed curve) and 
interactions out to next-nearest neighbors (NNN) (dotted 
curve). The transfer- matrix calculations were performed with 
5 steps, (L)=5, "(/(s) given by Eq. and fcBr = 0.84e, where 
e is the energy of a single kink. It appears that the nearest- 
neighbor-only interactions are in better agreement with the 
generalized Wigner distribution, although this judgment is 
hampered by the small system size and the breakdown, at 
large A, of the continuum step model. The dashed and dot- 
ted lines indicate limits of a^g for NN and NNN interactions, 
respectively, as derived from Gaussian approximations of the 
the continuum step model in the limit of large ^0,3. 



Interestingly, when all steps interact (via the potential 
given by Eq. the variance of a TWD with a fixed 
value of A will be smaller than the variance of a TWD 
with the same value of A and just nearest-neighbor steps 
interacting. Numerical evidence supports this intuition, 
and in fact it appears that the Wigner distribution de- 
scribes TWDs from nearest-neighbor interactions better 
than TWDs from all steps interacting (see Fig.O). On the 
other hand, the phenomenological approach which leads 
to the generalized Wigner distribution depends only on 
the interaction between nearest neighbor steps, which is 
identical in both cases. 

The difference appears to come from terms in Eq. ()25f) 
that are neglected in Eq. Since the variance is de- 

creased by interactions with more steps, it is clear that 
increasing the number of interacting steps will not in- 
crease the linear term V in Eq. H25II — that would in- 
crease the variance. Instead, it appears that increasing 
the number of interacting steps increases the quadratic 
terms in Eqs. (|19|l and H21|) . This would decrease the 
variance and correspond to the intuitive notion that the 
steps "interact more strongly" while still observing the 
constraint (s) = 1. Unfortunately, it also destroys the 
separability of the Schrodinger equation. 

A variety of circumstances, such as the presence of ad- 
sorbates |33| or electronic surface states [SJ, |33 , can give 
rise to short-ranged attractive forces between steps, os- 
cillating step-step interaction potentials, and other com- 
plicated interactions. Under such circumstances, one 
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FIG. 4: [Color online] A transfer-matrix calculation (•) of the 
TWD for a system of five steps with a mean step separation 
of 5 lattice units and the oscillating step-step potential given 
by Eq. 13711 . The curve is the solution of the corresponding 
Schrodinger equation and contains no free parameters. 



typically finds terraces of two more-or-less well-defined 
widths; interactions with steps beyond nearest neighbors 
can segregate these widths, leading to "step bunching" 
|36j |. This situation is somewhat analogous to liquid- 
gas coexistence, and as with liquid-gas coexistence the 
compressibility diverges. As a result, U{s) is linear in s 
[Eq. 

Vicinal surfaces with both elastic and surface-state me- 
diated electronic interactions may be characterized by 
potentials of the form |3 [s^ 

V{L) = AL-^ + BL-^''^ cos(2fcFi + (t>) , (36) 

where A is determined by the elastic interactions, B is 
determined by the couphng to the surface state, /cf is the 
Fermi energy, and is a phase shift. The oscillations of 
this potential can give rise to a coexistence as described 
above, together with a diverging compressibility and a 
linear U{L). As a dramatic — albeit rather unphysical 
— example, in Fig.0|we compare transfer-matrix TWDs 
derived from the corrugated interaction potential 

y(s) = 4cos(27rs) (37) 

[a special case of Eq. (|36ll ] with corresponding solutions 
of Eq. 0. Since the TWDs now have multiple peaks, 
clearly neither the generalized Wigner distribution nor a 
Gaussian distribution is appropriate. Again, considering 
that there are no free parameters, the phenomenological 
TWD is in excellent agreement with the transfer-matrix 
TWD. The slight disagreement is probably due to finite- 
size effects in the transfer-matrix calculation (which in- 
volved only five steps with {L) = 5). 



IV. SOLVING FOR THE STEP-STEP 
INTERACTION 

The chief practical interest in terrace width distribu- 
tions comes from a desire to use experimental data to 
determine the interaction potential V{L), or equivalently 
V{s). In this section we discuss three progressively more 
sophisticated methods. The first two are appealing in 
their simplicity, but they suffer from serious flaws, and 
therefore should be avoided. By building on them, we 
introduce the third method, which is more computation- 
ally intensive than the preceding two, but which has the 
flexibility necessary to be appplied to real experimental 
data. For more details see Ref. i37i 

A. Naive direct numerical approach 

In Sec. m we derived a Schrodinger equation that in- 
cluded V{s) in an obvious way, and in Sec. IIIII we have 
seen that similar Schrodinger equations yield TWDs for 
more complicated interactions that are in excellent agree- 
ment with numerical simulations of the TSK model. It is 
tempting to follow an analogous path when dealing with 
experimental data, using the following procedure: 

1. Connect the experimental data points with a 
smooth interpolating function P{s), such as a cubic 
spline. 

2. Calculate V(s) = ^/P{s). 

3. Find d^ij/ds^. 

4. The total potential is V{s) + uls) - E = 
[V'(s)]-Md'V'(s)/ds2]. 

5. At large s, V{s) is negligible, so the total is given 
by U{s) — E, which should be constant plus either 
a quadratic function of s or a linear function of s. 

6. Subtract U{s) — E from the potential to recover the 
step-step interaction V{s). 

This method has two serious problems. Experimental 
TWDs will contain measurements indicating P{s) = 0, 
which lead to division by zero in determining the total 
potential. Furthermore, typical statistical fluctuations 
make the numerical estimates of d^ip/ds^ from experi- 
mental data extremely unreliable. 

B. Fitting the TWD to a preconceived form of P{s) 

One may eliminate both of these problems by first per- 
forming a least-squares fit of the TWD to some positive 
definite and twice-differentiable function. For example, 
one might use 

P(s) = 005*^(1 + cis -I- C2S^) exp(-6s2) , (38) 



6 



where cq is a normalization constant and Ci, C2, and b 
all are parameters to be fitted. The fitted P{s) can then 
be analyzed as above. 

The shortcomings of this second approach are more 
subtle. Eq. (|38|l . for example, will always yield a poten- 
tial that diverges as as s approaches zero. In fact, 
we have fitted Eq. (|38|l to the Monte Carlo data shown 
in Fig. 121 Instead of correctly reproducing the potential 
V{s) = A3S~^, the analysis of the fitted P{s) produced an 
interaction that is well approximated by V{s) ^25"^. 
The intersection of this derived potential with the known, 
true potential is for a value of L near (L), which is not 
surprising — the TWD contains more information about 
the interaction where P{s) is large. 

C. Recommended: Fitting tlie interaction V{L) to 
a preconceived form 



0.15 




FIG. 5: [Color online] Monte Carlo (o) and phenomenologi- 
cal (solid curve) TWDs for steps interacting with a potential 
V'(s) = 2s~^ The dashed curve shows a least squares fit using 
a potential of the form given in Eq. II4UII . 



The method we do recommend is more computation- 
ally demanding, but it suffers from none of the above 
defects: 

1. Parameterize V(L), and make a crude initial esti- 
mate of the values of the parameters. At least some 
information on the form of the interaction is often 
available — whether it is influenced by a surface 
electronic state, for example. In principle, U{L) 
contains a term linear in L and a term quadratic in 
L, as in Eq. in practice, one can often use ei- 
ther a purely quadratic form of U{L) (as was done 
to derive the generalized Wigner distribution) or a 
purely linear form of U{L) (if the compressibility 
diverges) . 

2. Solve Eq. numerically for the specified param- 
eters to find ^'o(-Zj). 

3. Find 

x'^J2['^oiL)-Po.ALf (39) 

L 

for this set of parameters, where P(.xp{L) is the ex- 
perimental TWD. 

4. Iterate procedural steps 2 and 3 in a minimization 
routine to find parameters for which ^'o(i) best fits 
the experimental TWD. 




-20 ' ' 
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FIG. 6: [Color online] The solid curve shows the step-step 
interaction actually used (solid curve) to generate the Monte 
Carlo TWD shown in Fig. El The dashed curve shows the 
step-step interaction of the form given by Eq. 1401 with the 
parameters given in TableQ which were determined by a least- 
squares fit of the phenomenological TWD to the Monte Carlo 
TWD. To facilitate a better comparison, the fitted potential 
has been shifted by a constant energy so that it coincides with 
the true potential at s = l. 

In Fig. El we apply this approach to the Monte Carlo 
TWD shown in Fig.O We parameterize V{s) by 

V{s) = iis-i + A2S-^ + A^s-^ + A^s-'^ ; (40) 



TABLE I: Parameters in Eq. (I40II determined by a least- 



squares fit of the phenomenological TWD to the Monte Carlo 


TWD (see Fig. 














Ai 


A2 


Ai 








true 








2 









fit 


6.862 


1.1035 


0.40164 


0.05215 





this form amounts to a generalization of the two- 
parameter fit discussed in Ref. 138. where the only al- 
lowed term in V{s) was A2S~'^. The values of the fitted 
parameters are given in Table they clearly are quite 
different from those of the true potential. The resulting 
V{s) is shown in Fig. j^l where it is compared with the 
actual potential used to generate the Monte Carlo data. 
What is important, however, is that the fitted and true 
potentials are close over the range in s corresponding to 
significantly nonzero P{s) — in this case, 0.5 < s < 1.5. 
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The phenomenological TWD (Fig. [Sj) derived from these 
fitted parameters is in good agreement with the Monte 
Carlo TWD and is very difhcult to distinguish from the 
TWD derived from the true potential via Eq. H33|) . This 
result underscores the sensitivity of the process of ex- 
tracting V{s) and the necessity of very good statistics 
in experimental TWDs [ssj and a good idea of the func- 
tional form of V{s). 

For this reason, we strongly recommend making a si- 
multaneous fit of V{L) to experimental TWDs from sev- 
eral different misorientations, and preferably from sev- 
eral different temperatures. We are currently applying 
this approach to Cu surfaces vicinal to the (001) plane 
■ We have selected 10 experimental TWDs to be fit- 
ted, with mean step separations ranging from 5.44a to 
9.53a, where the lattice constant a = 0.255 nm, and with 
temperatures ranging from 285 K to 360 K. Our prelimi- 
nary studies |39l | show that even with experimental data 
of good quality, the extraction of real potentials is not a 
trivial matter. 

Using all 10 experimental TWDs, we were able to ob- 
tain good agreement between fitted TWDs and experi- 
mental TWDs; two examples are shown in Fig. [7| The 
corresponding interaction potential is shown in Fig. |S1 
clearly, it deviates markedly from Eq. The irregular- 
ities in V{L) make it plausible that a coexistence exists 
between two step widths, which would result in a diverg- 
ing compressibility. This is supported by the fact that 
the experimental TWDs are fitted with U (L) oc L than 
with U{L) oc L^; furthermore, the linear form of U(L) 
would lead to "fatter" tails, as were observed in Ref. 129. 

Although these results show internal consistency, the 
fitted interaction potential is markedly different from po- 
tentials derived from simulations of surface relaxation us- 
ing an empirical many-body potential 0i ^1 l4^ |43| . 
which are in reasonable agreement with Eq. We be- 
lieve that our interaction potentials are "over-fitted" to 
the experimental data, with the irregularities of the V{L) 
being too strongly influenced by statistical noise in our 
experimental data. This is not really surprising, because 
V{L) depends quite sensitively on P{L). 

A similar problem '3, '45] occurs in determining em- 
bedding functions for the Embedded Atom Model from 
experimental data; if too many parameters are fitted, the 
resulting embedding functions of course can fit the sup- 
plied experimental data, but they give worse agreement 
with experimental data not included in the fitting pro- 
cess. The recommended remedy 0, is to divide the 
available data into two groups, one of which is used to 
perform the fits, and the other of which is used to test 
the robustness of the fits. 

Presumably, the same remedy will work in our case. 
However, a more complete treatment of the problem, in- 
cluding estimates of the uncertainties of fitted functions, 
is deferred to a later paper [3^ . 




FIG. 7: Two experimental TWDs for Cu surfaces vicinal to 
the (001) plane. The circles are experimental data courtesy 
of M. Giesen, and the diamonds are experimental data cour- 
tesy of R. van Gastel. Both experiments were at 295 K. The 
corresponding fits were made by solving Eq. 1331 with the 
potential shown in Fig.|H]and U{L)(xL. 




FIG. 8: The step-step interaction potential for for Cu sur- 
faces vicinal to the (001) plane, as determined from fits to 
ten experimental TWDs like those shown in Fig. |7| This po- 
tential is in good agreement with experimental TWDs, but it 
has probably been "over-fitted." 



V. CONCLUSIONS 

We have seen that the generalized Wigner distribution 
can be derived from a Schrodinger equation somewhat 
in the spirit of the Gruber-MuUins approximation, and 
that straightforward extensions of this method work for 
general step-step interactions. 

In Sec. |n] we saw that the generalized Wigner dis- 
tribution can be derived exactly from the two-fermion 
Calogero model This, in turn, can be justified from 
a phenomenological model in which the force on two ad- 
jacent fermions is derived from the position-dependent 
pressure exerted by other fermions confined in a box, 
the size of which is presumably related to the correlation 
length in the x-direction. It is worth noting that the "en- 
tropic repulsion" is handled implicitly by the uncertainty 
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principle in the quantum mapping, so we need explicitly 
consider only the energetic interactions. 

Nothing in this phenomenological picture requires the 
step-step interaction to be given by Eq. ||SJ), and in 
Sec. mil wc demonstrate numerally that for very gen- 
eral step-step interactions the phenomenological picture 
yields TWDs in excellent agreement with numerical sim- 
ulations of the TSK model. This success is particularly 
impressive when a coexistence between two well-defined 
step widths occurs, as can happen when steps bunch 
[33 ■ Under this circumstance the compressibility di- 
verges, which causes the tail of the TWD to decay ex- 
ponentially with L rather than exponentially with , as 
it does in the case of the generalized Wigner distribu- 
tion. This phenomenon could explain the curious slowly 
decaying tails of TWDs mentioned in Ref. 25. 

Since exact solutions are available for both equations 
in the important case in which V{s) is given by Eq. ©, 
it is possible to find not only the TWD but also, using 
the methods of Ref. [s^, an improved estimate of step 
wandering. This will be undertaken in a separate paper. 

Finally, although we have addressed this work primar- 



ily to its surface science applications, it may be of in- 
terest to research in random matrix theory as well. The 
Calogero-Sutherland model with A — —1/4, 0, or 2 corre- 
sponds to random matrices with specific symmetries, and 
attempts to interpolate between them have simply varied 
the fraction of matrices belonging to each symmetry in 
the ensemble of random matrices. It is far from clear, 
however, what relation the Calogero-Sutherland model 
with more general values of A or our treatment of it has 
to such mixed ensembles. 
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